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Abstract 

There exists a well-known relation between the zeros of sine function, 
Bernoulli numbers and the Riemann Zeta function. In the present paper, 
we find a similar relation for zeros of q-sine function. We introduce a new 
q-extension of the Bernoulli numbers with generating function written in 
terms of both Jackson's q-exponential functions. By q-generalized mul- 
tiple product Leibnitz rule and the q-analogue of logarithmic derivative 
we established exact relations between zeros of sin, x and our q-Bernoulli 
numbers. These relations could be useful for analyzing approximate and 
asymptotic formulas for the zeros and solving BVP for g-Sturm-Liouville 
problems. 



1 Introduction 

One of the most impressive applications of Bernoulli numbers is related with 
zeros of sin a; function and the Riemann Zeta function. In the present paper, by 
proper generalization of Bernoulli numbers to g-BcrnouUi numbers, we estab- 
lish similar relation for zeros of sin^ x function. The generating functions of our 
Bernoulli polynomials and Bernoulli numbers are defined in terms of Jackson's 
q-exponential functions. Our generating functions and Bernoulli numbers are 
different from the known in literature [5], [5], [J. By q-differentiation of the 
generating function we get the recursion formula for q-Bernoulli polynomials 
which is reducible to the standard recursion formula in the limit g — > 1. From 
the definition of q-exponential functions we write the power series expansion of 
this generating function and obtain first few q-Bernoulli numbers. In order to 
write siuq x function as an infinite product in terms of its zeros, we introduce 
the q-generalized multiple product Leibnitz rule and the q-analogue of loga- 
rithmic derivative. This gives us relation between zeros of q-sine function and 
q-Bernoulli numbers. In the limit q ^ 1, our results transforms to the known 
relations between the zeros of sine function, Bernoulli numbers and Riemann 
Zeta function. 
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2 Zeros of Sine Function and Riemann Zeta Func- 
tion 



First we briefly review the known relation between the zeros of sin x function, 
Bernoulli numbers and the Riemann Zeta function. 

2.1 Bernoulli Polynomials and Numbers 

The generating function 

F,{z) = ^^ = f2B.^{x)^ (1) 

n—O 

determines the Bernoulli polynomials in x,Bn{x), \fn > 0. 

By differentiating in x we get the recursion formula for Bernoulli polynomials 

=nS„_i(x), n>l. (2) 

In addition, we have 

yn>l, Bn{x + 1) - Br,{x) = nx"-'^ . (3) 

Bernoulli numbers are defined as Bn{0) — bn- Then the generating func- 
tion for Bernoulli numbers follows from ([I} 



ri=0 

Below we display first few Bernoulli polynomials and numbers 

1 1 3 1 

Bo{x)^l, Bi{x)^x--, B2{x) ^ x'^ - X + B3{x) = x^ - -x'^ + -X. 

6o = 1, &i = ^2 = ^, 63=0. 

Bernoulli numbers are related with zeros of sinz function, and allows one to 
calculate the values of the Riemann Zeta function at even numbers argument 
[4] . We consider infinite product representation for sin z : 



00 

sin z = z 

n=l 



n i-iil p) 



cos 2 1 



sinz z ^ 1 



n— 1 Tz^n^ 



2 



^ n^-K^ 1 fij ^ n^TT 



n=l 



oo oo 2fc 

zcotz = l-2^^-^ (6) 

n=l fc=l 

From another side, we can represent this sum in terms of the BernouUi numbers. 
In the generating function Q 



n=0 

where b2n+i = for n, > 1, by choosing x = 2iz and 

2iz ze~*^ z(cosz — i sinz) , (2iz)" i, (2*^: 



g2iz _ I sinz sinz -^-^ n! (2fc)! 

n=0 k=l ^ ' 

we get 

zcotz = l-5]fe2fe(-l)^-^^. (8) 

Here we used the fact that b2k+i = for fc = 1, 2, .... It foUows obviously from 
observation that l.h.s. is even function of z. 

In this form, function on the Lh.s has infinite set of simple poles at z = 
±7r, ±2tt, .... If |z| < TT, then it is analytic and has unique expansion to Taylor 
series around z — 0. 

Comparing the expressions ^ and ([8]), we obtain 

°o 1 92fe-l 

This gives relation between different power of zeros of sine function and Bernoulli 
numbers. The left-hand side of this equation is the Riemann Zeta function 



oo 

c(^) -T.^ (10) 



n—l 

of even argument. This is why, we get expression of this Zeta function in terms 
of Bernoulli numbers 

C(2fc) = (-l)'=-i52.^7r2'=. (11) 
The following are the first few values of the Riemann zeta function: 

°° 1 2 
n—l 
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^(4) = y i- = — > (13) 

^ ^ n4 90' ^ ' 

n—1 

C(6) = E;;f-^- (14) 



n=l 



3 q-BernouUi Numbers and Zeros of q-Sine Func- 
tion 

Now we are going to find the similar relation between zeros of the g-sine function 
and the g-BernouUi numbers. 

3.1 g-Bernoulli Polynomials and Numbers 

First we introduce the q-analogue of Bernoulli polynomials and Bernoulli num- 
bers. The generating function for g-Bernoulli polynomials is defined in terms of 
Jackson's q-exponential functions as follows 



zeq(xz) zeq{xz)eq{—^) 
^.(f)(e,(f)-e,(-f)) =e,(|)-e,(-§) 



Y^BUx)^, (15) 



where the Jackson's g-exponential functions are [T] 



\n\„\ — ' \n 



and [n]q\ = [l]q[2]q...[n]g, [n]q — Two g-exponential functions are related 

to by the next formula 

eq(x)Eq{~x) = 1. 

By g-differentiation the generating function with respect to x, it is easy to obtain 
the recursion formula 

D-^Bl{x) = [n]qBl^^{x), (17) 

where Dqf{x) = and Bq{x) = 1. In the limiting g — 1, this relation 

reduces to the standard recursion formula ([2]). 

For n > 0, 6^ = i?^(0) we called the g-Bernoulli numbers. 

According to above definition, the generating function for g-Bernoulli num- 
bers is given by 



n 



4 



By the definition of Jackson g-exponential functions (jl6p we expand this 
generating function as 



i^.(f)(e,(f)-e,(-i)) 



[2]!+ 4 [4]! 



1 



+ 2 + (HsIT ^ + (2^W + 2%]t) + ^'^ (2* [2]! [3]! + 2%]! + FTsll) + 



1 



1 + f + Az2 + Bz3 + Cz4 . 



l-(^i + Az'+ Bz' + Cz' 



^+Az'' + Bz-" + Cz'^ + - + ^^'^ + ■B-z'' + Cz" 



+Az2^ 5^3 + Cz^ + ...)' + ..., (19) 



where 

A W 



c ^' 



22 [3]!' 
-7^ + 1 



23 [3]! ' 
_ [5]g6 + 1 g 



24 [5]! 24 [2]! [3]! 

Comparing terms with the same power of z we get first few g-Bernouhi 
numbers. 

For term we have 

4 2 [2]! 2 4 J y 

and 

For term z^ we get 

- E ^31- - [21^ + - X _ + 1 ) (20) 
^4- 24 \^^\- + [3], [2]! [5] [4] ) ^^"^ 

and as a resuh 

4" 24 l^^^J- + [3]! [2]! [5] [4] 
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i.S = l, l>t = --. l,l = -l^2]-^^-gj, 6? = 0, (21) 



2^ y ' ' [3]! [2]! [5][4] 
By choosing z = 2it in generating function (jlSp . we obtain 



Eq{it) {eq{it) - eq{-it)) Eq{it)smqt siugt 
From the q-analogue of Euler identity eqlix) — coSq x + i sin^ x, we have 

^o(2ji) = (cos„ t — i sin„ = t cot„ t — it = > 

sing i In I' 



n=0 



6^i + 6?(2zt) + f;6^^^, 



where 

cos„ t 

cot„ t = — 

sinq t 

and 

_ eg(it) + eq{~it) . eq(it) - eq{-it) 

' 2 * 2i 

Then, substituting and 6^' into the above equahty we get 



tcotqt=l + f2bl^^, (24) 



tcot,t = l + ^6^,i||^. (25) 

k=l ^ ^' 

Here the left-hand side is even function of t, so that in the last sum odd coeffi- 
cients vanish 62fc+i = for fc = 1, 2, .... 



3.2 Zeros of g-Sine Function 

Now we like to express the l.h.s. of ([^5]) in terms of zeros of sin^ x function. We 
start with proposition : 

Proposition 3.2.1 q- Generalized Multiple Product Leibnitz Rule: 

Dq{h{x)f2{x)...fn{x)) - {Dqhix)) h{x)...U{x) 

+ h{qx){Dqh{x))h{x)...U{x) 
+ ... 

+ h{qx)f2{qx)...fn-i{qx){DqUx)) (26) 
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Proof 3.2.2 For n ~ 1 it is evident. For n = 2 it gives q-Leibnitz rule |7]/ 

Dg {h{x)f2{x)) = {D^h{x)) ,Mx) + hiqx) {D^hix)) . 
Suppose it is true for some n. Then by induction 

Dq{h{x)h{x)...fn{x)fn+l{x)) = D ,{h{x) J2{x) .. . Ux)) fn+l{x) 

+ fl{qx)f2{qx)...fn{qx) {DqJn+l{x)) 

= {{D,h{x)) ...Ux) + ... + h{qx)... {D,fn{x))) fn+i{x) 

+ fl{qx)f2{qx)...fn{qx) {Dqfr,+ l{x)) , 

which is the desired result. M 

According to the above proposition we have the foUowing rule of differentiation 
multiple product (the g-analogue of logarithmic derivative) 

D,{flf2...fn) ^ mx) h{qx)f^{x) h{qx) fn-l{qx) f'^(x) 

hh-fn h{x) ^ h{x) /2(a;) h{x) - /„_i(a;) Ux) ^ ' 

Example: If fk — (x — Xk) and fi...fn = Y[k=i(.^~^k) is function with n zeros, 
xi,...,Xn, then we have 

Dq([\k=ii^ ^ ^k)) _ 1 , (qx-xi) 1 [qx - xi) {qx - X2) 1 



]Xk=i{^-^k) {x - xi) {x - xi) {x - X2) [x - xi) {x - X2) (x - zs) 

_^ ^ [qx - xi) {qx - X2) {qx - Xz) {qx - Xn-i) 1 

{x - Xl) {x - X2) {x - X3,) {x - Xn-l) {x - Xn)' 

as a simple pole expansion. 

Expanded to simple fractions this expression can be rewritten as 



nLi(2^-a;fe) tr[x-xk 



E^^' (28) 



where coefficients 

Dq{Uk=li^ - Xk)) 



Ak = Resi^^- 
= Resi^^ 



1 {qx — Xl) 1 {qx — Xl) {qx — X2) 



{x - Xl) {x - Xl) {x - X2) {x - Xl) {x ~ X2) {x ~ Xn) 

Particularly, for n = 2, 

1 qx — Xl 



Ai = lim (x - Xl) , , 

x^xi \ \X — Xl {X — Xi){x — X2) 

= l + = (29) 

{Xl — X2) Xl - X2 
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A2 = lim (x - X2) 



1 qx — Xi 



x->-X2 \ \X ~ Xl (x Xl){x — X2) 

_ qx2 - Xl 

X2 — Xl 

and we get 

Dq{{x ~ Xi){x ~ X2)) _ f qxi ~ X2\ 1 f qx2 - Xl 



(30) 



{x — Xi){x — X2) \ Xl — X2 J X — Xl \ X2 — Xl J X — X2 

We consider sin^ x function as an infinite product in terms of its zeros 
Xn = Xn{q) in the following form 



(31) 



sin,x = xl[( I- -^\=x( I- ■-,)(!- ■-,)... (32) 

n— 1 ^ ^ ^ ^ 

By using the above property ([?f|. we have 



D„sin„a; ^9 11^^=1(1 



siuo X 



COtg X ~ 



1 a.(i-M(-M 



I- 



,33, 



X (^^ 



where we ordered zeros as |a;| < \xi\ < \x2\ < ... < |a;„| < so that |^| < 1, 
for any k. The above expression can be written in a compact form as follows 

00 n-1 M _ (j2 2^ 

X cot, x^l-[2]qY: jf^ n V (34) 
3.2.1 Quadratic order: 

Now we compare expressions psp and p4p by equating equal powers in : 



2 [2]! 4 [4]! 
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1 + (1 - g^) + (1 _ 



2\ 2 



^2 



2 / 2 \ 2 

2\ I /I 2\ ' 



l + + +••• • (35) 

At the order we have 



n=l 



and using (PT|) for the value of Bernouhi number 6| = i ^[2] — g — -jij^ , we 
obtain 

°° 1 1 

In the hmiting case q — > 1, [3]! = 6 and we have 

°° 1 1 
hm y = — . 

Due to relation 



4c(2) = E 

it implies 



1 1 



n^TT^ 6' 

n— 1 



lim Xn{q) = niT. 

We found relation between zeros Xn of g-sine function and 63 at order x'^. 
3.2.2 Quartic order: 

Now we will find relation at the order x*, this why let us call 

^2 



2 — 



then the above expression is written in terms of ^ as follows 

X cot, - = 1 - [2]g E n = ^ ~ ^^4t"' + ''m"' + - (''^ 

For simplicity we denote 



n=l ' fe=l 
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then open form of the above expression gives 

1-6 1-6 1-6 1-6 1-6 1-6 
, , 6 (1 - g^6) (i-g^6-i) , 

1-6 1-6 1 - 6-1 
For = 161 < 1, Taylor expansion of the above expression is 

A = 6(i + 6+?? + -) + 6(i + 6 + 6' + -)(i + (i-'?')6 + (i-'7')^? + -) 

+ 6(1 + 6 + CI + •••) (1 + (1 - g')6 + (1 - q^)ei + •••) (1 + (1 - 'z')6 + (i - a')6' + •••) 
+ ... + 6(1 + 6 + ?^ + ...)...(i + (i-a')6-i + (i-a')CLi + •••) + ••• (39) 

Here we should consider just terms to collect order x'^, so we denote 

B = + 6' + ••• + en + ••• + 66(l-9') + 66(l-«') + 66(l-9') + ••• 
+ 66(1 - + 66(1 - + ... + 66-i(i - + - 

oo 

= J2ek + {l-q')C, (40) 
fe=i 

where 

oo oo oo 

C7^^66 + E^26 + -+ E + - (41) 

fc=2 fe=3 A;=n+1 

By 

n 

6 = 'S'n 

fc=l 

and 

lim Sn = S, 

n—>oo 

then we can write the sums as 

oo 

k=l 

oo OO 

A;=2 k=l 
oo oo 

fc^3 fe=l 
oo 

E^fe = ^-^"-i (42) 



fe=n 
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Rewriting (PO]) in terms of S, we obtain 

B = S^ + {l-q^){^,{S-S,) + ^2{S-S2) + ... + US~S,,) + ...) 

- 52 + (1 - _ ^^^^ _ ^^(^^ ^ _ _ ^^(^^ ^^^^ ^ 

= ^' + (l-9')(5'-E^fe + ^), (43) 



where 



or 



^ = -6^1 - + 6) - - - + 6 + - + ^n-l) + ... (44) 



D = -6(6+6 + ...+en)-6(e3+e4 + ...+e«)-... -^«(en+l + ...)-... 

= -US-Si)-US-S2)-...-US-Sn). (45) 
Comparing with (PT|) we find D ~ —C, then by equating (HUl) and 

5 = ^2^(1-^2) Is^-f^ek-c] =S' + {l-q')C, (46) 



fc=i 



we get 



It gives 



1 1 °° 

C-ls'-lY.^l (47) 



2 2 

fc=i 



00 / ^ ^ 00 N 

- E^' + (i-'^')U^'-^E^' 



fe=i 



For term then we have 

bl^x' = -[2]qB, (49) 
and substituting — ^ and S = J^kLi = X]fc°=i finally we obtain 



k 
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where 



4 2« V'"" ' [3]! 12]! |5||4| 

In the limitmg case — > 1, 



°° 1 1 
Km > — -— — = — . 



90- 

From the relation (fT3)) we get 

1 °° 1 1 

n— 1 

This procedure can be continued to higher order relating higher powers of 
zeros of q-sine function with q-BernouUi numbers. However, procedure at every 
next step become more and more complicated. Moreover, it is another problem 
finding exact form of zeros x„(q) of siug x from infinite set of equations. 



4 Conclusions 

In this paper, we have introduced the set of q-BernouUi numbers and found 
relations between these numbers and zeros of g-sine function. The exact formula 
for these zeros is not known, this is why, our relations can be useful in analyzing 
asymptotic formulas and approximate form of the zeros. 

Here we mention some simple approximate formula for zeros of q-sine func- 
tion. It is coming from analyzing graph of siug x and for g > 1 is 



then, 



X2=q^xi, xs. ^ q'^X2 ^ q'^Xi,...,Xn = q^^" ^^Xi, 



1 1 1 1 1 



[3]! xl x\ xl ■ x\ q^x\ q^x\ ■ 



and from sum of geometric series in ^ we have 

1 1 (?4 



[2] [3] x\q^ 
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From this expression we have the first root as 



As a result, ([5^ can be written in the following form 



sm^ X = X 

n=l 



where for wave number we have the discrete set — pj^^rq-- 

These results can be useful in solving BVP for q-wave equation and g-Sturm- 
Liouville problem. These questions are under investigation now. 
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